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Abstract. The plan for mass and spin measurement of SUSY particles with the ATLAS detector 
is presented. The measurements of kinematical distributions, such as edges in the invariant mass 
of leptons and jets, could be used to constrain the model of SUSY that may be discovered at the 
LHC. Examples from a few points in the mSUGRA scenario are provided with an emphasis on 
measurements that can be conducted within the first few years of data taking. 

PACS. ll.30.Pb Supersymmetry - 12.60. Jv Supersymmetric models 



1 Introduction 

Discovering Supersymmetry (SUSY) is one of the mo- 
tivations for building the Large Hadron Collider (LHC) 
which is scheduled to take data in 2008. ATLAS (A 
Torodial LHC ApparatuS) is one of the two general 
purpose experiments at the LHC. Currently the AT- 
LAS experiment is revisiting its SUSY studies in the 
framework of the CSC ("Computing System Commis- 
sioning"). For this exercise event samples for a set of 
SUSY benchmark points were produced with the de- 
tailed simulation of the ATLAS detector. The aim of 
the work is to better understand the impact of realis- 
tic experimental conditions including trigger efficien- 
cies and imperfect calibration/alignment on the SUSY 
potential of the experiment. 

If SUSY is discovered by inclusive searches, next 
step will be to measure the masses and the spins of 
the SUSY particles through the analysis of exclusive 
decay channels. In this note SUSY mass and spin mea- 
surement techniques and methods for determining un- 
derlying SUSY model parameters are described within 
the mSUGRA framework. 



2 mSUGRA Framework 

The minimal SUSY extension of the SM (MSSM) brings 
105 additional free parameters into the theory thus 
making a systematic study of the full parameter space 
difficult. A specific well-motivated model framework 
is usually assumed in which generic signatures can 
be studied. In the minimal Supergravity (mSUGRA) 
framework, SUSY is broken by the gravitational inter- 
actions and the masses and couplings are unified at 
the grand unified energy scale giving five free parame- 
ters; mo and mi/2 (universial scalar and gaugino mass 
parameters), Aq (universial trilinear coupHng), tan(/3) 
(the ratio of the vacuum expectation values of the 
scalar fields) and sgn(/i) (the sign of the higgsino mass 



term). In the R-parity conserving mSUGRA models, 
SUSY particles are produced in pairs and the lightest 
SUSY particle (LSP) is stable and undetected result- 
ing in large missing transverse energy (dominant sig- 
nature). It will be possible to detect SUSY in multiple 
inclusive signatures over a large part of the mSUGRA 
parameters space with just 1 fb~^ of integrated lumi- 
nosity, as discussed in [Ij. As a part of the ongoing CSC 
exercise, several mSUGRA points that are favored by 
the WMAP data [2j have been simulated with a de- 
tailed detector simulation. The mSUGRA parameters 
for these points are given in Table [H 

3 Mass Measurement Techniques 

In the R-parity conserving mSUGRA models all SUSY 
events contain two invisible neutralinos (x?) which, be- 
ing only weakly interacting, escape the detector; there- 
fore no mass peaks can be reconstructed directly. How- 
ever kinematic endpoints and thresholds in the invari- 
ant mass distributions of the visible decay products 
can be measured. The values of these kinematic fea- 
tures can be expressed as a function of the masses 
of the involved sparticles [3j. The mass measurement 
strategy is to exploit kinematics of long decay chains 
originating from gluino or squark production. Specifi- 
cally, from the reconstruction of the decay chain ql 
QX2 ~^ Q^R^^ ~^ QXi^^l^ can measure the masses 
of Ir^ X2 Xi- Examples of mass measurements 
using the endpoint method are given below. 

3.1 Dilepton Endpoint 

A clear endpoint is exhibited by the dilepton invari- 
ant mass distribution for the decay channel X2 ~^ 
T^P Xil^P' Fully simulated samples for the SU4 
point (see Table [T]) and tt background are considered. 
Events are selected by requiring at least two leptons 
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Table 1. mSUGRA points chosen for full detector simulation in the CSC exercise. IS A JET 7.71 is used and the top 
quark mass is set to 175 GeV. 



Point Name 


mSUGRA Region 


mo (GeV) 


mi/2 (GeV) 


Ao (GeV) 


tan(/3) 


sgn(/i) 


a (pb) 


SUl 


Coannihilation 


70 


350 





10 


+ 


7.43 


SU2 


Focus 


3550 


300 





10 


+ 


4.86 


SU3 


Bulk 


100 


300 


-300 


6 


+ 


18.59 


SU4 


Low mass 


200 


160 


-400 


10 


+ 


262 


SU6 


Funnel 


320 


375 





50 


+ 


4.48 


SU8.1 


Coannihilation 


210 


360 





40 


+ 


6.44 


SU8.2 


Coannihilation 


215 


360 





40 


+ 


6.40 


SU8.3 


Coannihilation 


225 


360 





40 


+ 


6.32 



I Dilepton Edge invariant mass for SU4 

^ 20r \ \ \ \ r- 



■f 1 ndf 


22,57/ 29 


Prob 


0.7957 


pO 


49.17 + 2.876 


p1 


-0.2248 + 0.04664 
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Fig. 1. The dilepton invariant mass distribution for a full 
simulation sample of the SU4 point with an integrated lu- 
minosity of 0.35 fb~^. A triangular function convoluted 
with a Gaussian is fitted to estimate the edge position. 
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Fig. 2. The invariant mass distribution of lepton-jet for the 
SUl full simulation sample with an integrated luminosity 
of 5 fb~^. A triangular function convoluted with a Gaussian 
is fitted to estimate the edge position. 



(e or /i) with pt >10 GeV and |7^| <2.5, a calorimet- 
ric energy deposit Et <5 GeV in a 7^ — cone of 
size 0.3 around the lepton direction, and at least four 
jets with Et >100, 50, 50, 50 GeV and the missing 
transverse energy E^^^^ >120 GeV. The effective mass 
variable Meff, defined as Meff = E^'^^ -h J2}et^T,iet 
taking into account the first four leading jets, is re- 
quired to be Meff >550 GeV. The combinatorial back- 
ground from tt and decays cancel in the combina- 
tion e+e~ — e^/i^ (so called flavor subtraction) 
which is plotted in Figure [H The dilepton invariant 
mass distribution is shown for an integrated luminos- 
ity of 0.35 fb~^. The mass distribution is fitted to a 
triangular function convoluted with a Gaussian to esti- 
mate the edge position. The endpoint fit gives a value 
of (49.2±2.9) GeV which is within 1.6a of the expected 
value of 53.7 GeV. The signal significance is calculated 
as 16.5 for 100 pb~^ of data. 



with the lepton and then subtracting the mixed-event- 
jet distribution from the same-event-jets distribution 
(with a normalization correction applied) to obtain 
an inferred "correct jet" distribution. Fully simulated 
samples for the SUl point and tt background are con- 
sidered. Events are selected by requiring one lepton (e 
or /i) with Pt >20 GeV and |7^| <2.5, and a lepton 
isolation cut of Et <10 GeV in a cone of 0.45. Ad- 
ditional cuts are applied to reduce the tt background; 
the leading and second leading jets are selected with 
Et >200 GeV, the transverse mass Mt <60 GeV or 
Mt >100 GeV, E^'^^ >250 GeV. The invariant mass 
of the lepton-jet is shown in Figure [2] for an integrated 
luminosity of 5 fb~^. The mass distribution is fitted 
to a triangular function convoluted with a Gaussian 
to estimate the edge position. The endpoint fit gives 
a value of (283.6±4.8) GeV compared to the expected 
value of 284 GeV. The endpoint can be determined 
with 5(7 statistical signal significance for 5 fb~^ of data. 



3.2 Lepton-Jet Endpoint 



For the SUl point, the decay ql qxf ql^^i 
ql^h'iXi gives rise to an endpoint in the lepton-jet in- 
variant mass distribution. A mixed event technique 
is used to subtract the combinatorial jet background. 
The technique makes use of randomly pairing the jets 
from a different event (satisfying same event selection) 



3.3 Ditau Endpoint 

The mass difference between X2 Xi can also be 
measured from the endpoint of the ditau invariant 
mass distribution from the decay channel X2 -^t^t^ - 
XiT^r^. The SU3 point has a factor of four larger 
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X^/ndf 


8.507/ 12 


Prob 


0.7443 


pO 


2282 ±172.9 


P1 


3.993 + 


0.031 


p2 


0.2656 ± C 


.0234 
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Table 2. Determination of mSUGRA parameters for the 
SPSla point using the SFitter program. sgn(/i) is fixed. 
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Fig. 3. The ditau invariant mass distribution for the SU3 
fast simulation sample with an integrated luminosity of 10 
fb~^. It is fitted to a function as decribed in the text to 
estimate the endpoint. 



branching ratio to X2 ~^ Xi^^'^^ than to e//i. A fast 
simulation sample for the SU3 point together with 
Z+jets, W+jets, tf, 66+jets, dijets and multijets back- 
grounds is considered. Events are selected by requiring 
at least one jet with Et >220 GeV, at least three jets 
with Et >50 GeV and at least four jets with Et >40 
GeV and E^'^^ >230 GeV. The angular seperation 
between the taus is chosen as AR{r,r) <2. Only the 
hadronic tau decays are considered. The ditau invari- 



ant mass distribution calculated from r+r~ — r^r^ 
(to cancel the background from xt decays) is shown 
in Figure [3] for an integrated luminosity of 10 fb~^. 
Since the measurement of the endpoint by a linear fit 
has a strong dependence from the fitting range or bin- 
ning of the distribution, a new approach is considered 
to estimate the endpoint. A fit function in the form of 





SPSla 


^masses 


hedges 


mo (GeV) 


100 


3.9 


1.2 


(GeV) 


250 


1.7 


1.0 


tan(/3) 


10 


1.1 


0.9 


Ao (GeV) 


-100 


33 


20 



in [6]. The potential for extracting the masses from 
the endpoint measurements is described in [9j for the 
SPSla point [lOj with 300 fb-^ The ultimate goal is 
to perform a fit of the mSUGRA parameters from a 
given set of mass measurements using the tools such 
as Fittino [7J and SFitter [8j. As an example, the mea- 
sured masses and kinematical edges of the SPSla point 
from [9] are fed into the SFitter program to determine 
the mSUGRA parameters for 300 fb"^ [II]. As seen m 
Table[2]the parameters can be determined with a preci- 
sion at the percent level by using the masses {Amasses 
column), moreover the precision can be improved sig- 
nificantly by using the measured edges, thresholds and 
mass differences (Aedges column) in the fit instead of 
the masses. 



5 Spin Measurements 

If SUSY signals are observed at the LHC, it will be vi- 
tal to measure the spins of the new particles to demon- 
strate that they are indeed the predicted super-partners. 
Two methods of measuring the spin are described be- 
low. 



Po 

y = — .exp 

X 



2pl 



{ln{x) -pif 



(1) 5.1 Neutralino Spin Measurement 



is used (modified from [4J). The inflection point is cal- 
culated as: 



xip = exp 




(2) 



The measurement of the inflection point is more sta- 
ble under variations of the fitting range or binning. A 
calibration line is driven between the inflection point 
and the endpoint by varying the involved masses in 
the decay chain. This calibration line is found to be 
y = (0.47 ± 0.02)x H- (15 dz 2). The endpoint value is 
then calculated as (105 ± 4) GeV compared to the ex- 
pected value of 98.3 GeV. 



4 mSUGRA Parameters Determination 

More endpoint measurements can be done for instance 
from the reconstruction of gluino decays; g ^ X2^f 
and g ^ tit tbxf for which the details can be 
found in [5] and [6] respectively. Also, the mass of qr 
can be measured from the decay Qr ^ XiQ studied 



IT 

ar near 



The decay chain qX2 ^ Q^R^near ^ QXi^fc 

provides a good oppurtunity to measure the spin of X2 
using the lepton charge asymmetry [12 • The squarks 
and sleptons are spin-0 particles and their decays are 
spherically symmetric, X2 however has spin- 1/2 and 
the angular distribution of its decay products is not 
spherically symmetric. This leads to a charge asym- 
metry of the invariant mass of the quark and the near 
lepton which is defined as: 



dcF 



dm{ql^) 



(3) 



The asymmetry is suppressed by the fact that quark 
jets cannot be experimentally distinguished from anti- 
squark jets at the LHC, however at the LHC much 
more squark than antisquark will be produced, and 
therefore a residual asymmetry is still observable. This 
asymmetry is calculated for the SU3 point considering 
a fast simulation sample of 30 fb~^ together with the 
most relevant SM backgrounds, tf, W+jets, Z+jets. 
Events are selected by requiring two opposite sign lep- 
tons (e or fi) with pt >10 GeV and |7^| <2.5, a lep- 
ton isolation cut of Et <10 GeV in a cone of 0.2, 
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Fig. 4. The lepton charge asymmetry in the invariant mass 
of the lepton-jet for the SU3 point with an integrated lu- 
minosity of 30 fb~^. 
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Fig. 5. cosOii distribution for the Point 5 with an inte- 
grated luminosity of 200 fb~^. 



at least four jets with Et >100, 50, 50, 50 GeV and 
^miss ^200 GeV. The charge asymetry is plotted in 
Figure [H by considering both the near and far leptons 
since these are not distinguishable for the SU3 point. 
The confidence level for the asymmetry distribution 
to be flat (spin-0) is calculated using two independent 
statistical methods as described in [13j and found to 
be ~ 10~^. It is observed that 10 fb~^ of data would 
be sufficient to detect a non-zero charge asymmetry at 
the 99% confidence level for the SU3 point. 



5.2 Slepton Spin Measurement 

The spin of the slepton can be measured using an an- 
gular variable which is sensitive to the polar angle in 
direct slepton pair production qq Z^/j 
Xi^^Xi^" as discussed in The angular variable is 
defined as cosOii = cos{2tan~^e^^^i'^i- ^'^^) = tank 
{Ar]i+i-/2) It is interpreted as the cosine of the po- 
lar angle between each lepton and the beam axis in 
the longitudinally boosted frame in which the pseudo- 
rapidities of the leptons are equal and opposite. This 
variable is on average smaller for SUSY than for the 
Universial Extra Dimensions (UED) so it can be em- 
ployed as a spin-discriminant in slepton/Kaluza-Klein- 
lepton pair production in hadron colliders. The study 
was performed on Point 5 studied for the ATLAS TDR 
[15j which has a phenomenology very similar to the one 
of point SU3. A fast simulation sample of 200 fb~^ is 
considered together with the major SM backgrounds. 
Events are selected by requiring two opposite sign and 
same flavor electrons or muons with prih) >40 GeV, 
Prih) >30 GeV, Mu <150 GeV, no jet with Et >100 
GeV, no b-jets, E^'^^ >100 GeV, Mt2 <100 GeV. 
The angular variable is plotted in Figure [5] together 
with the predictions for SUSY (black line), phase space 
(dotted blue line) and the UED (dashed red line). The 
yellow shaded band represent the SUSY expectation 
when all the SUSY particle masses are simultaneously 
changed by ±20 GeV. As seen in the figure, the data 
points are much better matched to the slepton angu- 
lar distribution than to either the phase-space one or 
the UED-like one. This shows that cosOfj does indeed 



measure the spin of the sleptons for this point. Fur- 
ther studies with the Snowmass points SPSla, SPSlb, 
SPS3, SPS5 [lOj allow slepton spin determination with 
100-300 fb-i of data. 



6 Conclusions 

LHC brings experimental physics into a new territory. 
ATLAS can discover the SUSY particles up to 2-3 GeV 
mass scale if they exist. Understanding the detector re- 
sponse and the SM background will be a challenge at 
the beginning. Many techniques have been developed 
to measure the masses (edges, thresholds, mass differ- 
ences) and spin of SUSY particles and to determine 
the underlying model parameters. 
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